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(57) Abstract: Methods and apparatus for obtaining a repre- 
sentation of the distribution of electrical impedance within a 
multiphase flow with an electrically continuous or discontinu- 
ous principle flow (3) contained within an electrically conduc- 
tive solid ring electrode (1), comprising providing a plurality' 
of mutually spaced electrical contacts (2) mounted at the out- 
side wall of the ring and electrically contacted with the ring, 
applying currents or voltages (4) to the ring from the electrical 
contacts (2), generating a more homogeneous electric field dis- 
tribution within the material (3), measuring voltage or current 
(5) distribution along the ring from other electrical contacts 
(2), relatively intensifying the imaging sensitivity at the central 
area of the sensing domain using a 7t/2 angle sensing strategy 
and reconstructing the representation of the impedance distri- 
bution using CG method with an error processing strategy. 
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Electrical Impedance Tomography 

This invention relates to a method of measuring the distribution of electrical 
impedance of a multi-phase flow using an electrically conductive ring electrode, and 
5 to a method of detennining the internal structure of a body which utilises such a 
method and a computer programme product adapted for determining the internal 
structure of a body. 

The invention further relates to an apparatus adapted to measure the distribution of 
10 electrical impedance of a multi-phase flow using an electrically conductive ring 
electrode. 

This invention relates in general to electrical impedance tomography (EIT), and in 
particular to a new sensor construction and apparatus, a sensing strategy and an 
1 5 algorithm for measuring flow with electrically continuous and discontinuous 
conductivity phase, such as the oil/gas/water multiphase flow and the bubbling 
formation or foam construction. 

A number of sensor construction methods for EIT were proposed in previous patents, 
20 which were comprise of placing an array of spaced electrodes in direct contact with 
the human body or process fluid [GB2119520, GB2160323, US4617939], or in 
indirect contact with the human body coupled with electrolyte [GB20 19579], a 
sponger [WO98/23204] or any other suitable liquid medium [US448683 5] with 
conductivity matched to the typical expected conductivity of the object being 
25 measured. To obtain better distinction between a known conductivity and an 

unknown conductivity change in a body, a sensing method of the best pattern of 
currents was disclosed in US4920490 and US5272624, and a sensor arrangement 
using a number of electrically conductive compartment electrodes each separated 
from neighbouring compartments by an insulating partition was disclosed in 
30 WO98/23204. It has become known to apply EIT to vessels and pipelines made of 
electrically non-conductive material, such as acrylic or other plastic materials, or 
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made of electrically conductive metallic materials (WO 95/24155). la the process 
applications of EIT, it was also known that the electrodes were mounted inside the 
vessels or pipelines, and electrically insulated from one another, so as to be directly 
in- electrical contact with the contents within > 

5 

It is clear that previous EIT methods were based on utilising a number of discrete 
electrodes mounted at spaced locations of the wall, which had to be electrically 
insulated from one another and in electrical contact with body or material contained 
within the wall whether die wall is electrically non-conductive or conductive. The 1 
1 0 coupling materials used between the electrodes and the human body in the indirect 
contact method from previous knowledge, such as the- electrolyte and sponger, were 
purely purposed as an electrical coupling medium with a conductivity value 
comparable to that of the volume to be measured, which were neither suitable for * 
flow measurement nor functioned as a part of sensing strategy for improving the 
15 measurement capability. The limitations from the discrete electrode structure have 

been observed in following aspects. Firstly, it can not provide a correct measurement , 
from a fluid with a discontinuous phase of electrical conductivity since some of these 
electrodes have lost their electrical contact with the conductive phase of the fluid, e.g. 
a half filled horizontal stratified flow or a foam construction with a large bubbling 
20 formation [M. Wang and C. X Young, Chem. Eng. Science, 54(5), 1999]. Secondly, 
it generates a very non-uniform distribution of the imaging sensitivity because high 
dense currents are always present at the domains near to the current injecting and 
sinking electrodes. The limited contact area of the electrodes also produces a large 
common voltage, which highly reduces the signal-to-noise ratio for the differential 
25 measurement adopted in most cases, hi addition, it is difficult to present data with 
high repeatability from a sensor with discrete electrodes since the' unstable 
. performance of the contact impedance or the electrode-electrolyte interfaces. Finally., 
the intrusive construction of the electrodes in these sensors has obvious difficulties 
from its installation for satisfying most industrial environments. 
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The major features of this invention are to employ an electrically conductive ring 
with a conductivity value higher than that of the principle flow [Figure 1] as a 
tomographic sensor to supply a continuously electric current flow to the sensing 
domain, generate a more homogeneous imaging sensitivity distribution over the 
5 sensing domain and provide a smooth sensing surface without the obstruction effects 
to the flow. A PI/2 sensing strategy as an additional part of the invention further 
improves the imaging sensitivity distribution and an image reconstruction algorithm 
accompanied with the invention provides a practical solution for eliminat ing 
hardware noise in EIT applications. 

10 

Brief Description 

The object of the invention is to provide a sensing method and apparatus with an 
improved imaging sensitivity distribution for measuring the electrical impedance 
1 5 distribution in a principle flow with an electrically continuous or discontinuous phase 
and a practically applicable accurate reconstruction algorithm for flow measurement. 

Thus according to the invention we provide a method of measuring the distribution of 
electrical impedance of a multi-phase flow, of which the principal flow has an 
20 electrically continuous or discontinuous phase, using an electrically conductive ring 
electrode, characterised in that the ring electrode uses a PI/2 sensing strategy. 

The ring electrode may have a muchlugher conductivity-compared to-the : principle , 
flow or materials to be measured and it is employed as said a ring electrode of a part 
25 of the sensor body. A more homogeneous electric field distribution can be produced 
if the conductivity of the electrical conductive ring is much higher than that of the 
target content The electrically conductive ring may be made of solid substances 
such as a kind of metal, conductive rubber or ceramics. 

30 A number of electrical contacts may be embedded into the conductive ring electrode 
in good electrical contact with the outside wall of the ring electrode or penetrated 
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through the wall of the ring electrode in good electrical contact with the wall. The 
size of the contacts can be very small and the amount of the contacts can he installed 
to be much more than those in previous sensors. The sensitivity distribution is 
relatively intensified at a particulate interested area in the whole sensing domain, 
5 based on a specific sensing strategy or a combination of a basic set of independent 
measurements. 

The sensor geometry and its conductivity can be optimized for the objective that the 
condition number of the sensitivity matrix becomes less and the sensitivity field 
1 0 bocomes more homogenous. 

The method of the invention may be especially useful in determining the structure of; 
for example, process materials, such as a process fluid or a part of a human body, 
with electrically continuous phase such as that in a miscible liquid or discontinuous 
1 5 phase such as that in a foam formation or a half filled horizontal oil/water/gas 
stratified flow. 

The method may be extended to a sensor with discrete electrode structure using an 
external resistor network connecting all electrodes. Furthermore a multi-step image 
20 reconstruction algorithm using an error processing method is employed to reconstruct 
the conductivity distribution. 

According to a yet iurther feature of the invention we provide a computer programme 
product adapted for determining the internal structure of a body with electrical 
25 impedance tomography characterised in that the programme includes the use of the 
algorithm 

AVj ^ (to«°* J - 1 2 - P > (Equation 6) 

30 . ** 

Where/ is the measurement-projection location and k is the pixel number, s Jk denotes 
the sensitivity coefficient at pixel £ under the measurement-projection j, P denotes 
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the maximum number of measurements, w denotes the maximum number of pixels, 
a k and Acr k are the conductivity and conductivity change at pixel k 3 respectively, and 
Vj and AVj are in respect to "the reference voltage and the Voltage change at 
measurement-proj ectiony. 

5 

The solution of above equation can be implemented by assembling an inverse matrix 
using finite-step- CG method so that the on-line image reconstruction only involves 
the multiplication of a matrix and a vector. 

10 We also provide an apparatus adapted to measure the distribution of electrical 
impedance of a multi-phase flow of which the principal flow has an electrically 
continuous or discontinuous phase, using an electrically conductive ring electrode 
characterised in that the ring electrode is adapted to use a PI/2 sensing strategy. 

15 Detailed Description 

The invention will now be illustrated by way of example only and with reference to 
the accompanying drawings, in which 

20 Figure 1 illustrates the basic construction of the invention; 

Figure 2 illustrates a traditional EIT sensor under a stratified flow; 
Figure 3 illustrates a flange construction of the sensor; 

. Figure 4 shows two simulated electric field distributions with equi-potential lines 
from the sensor with discrete electrodes; 
15 Figure 5 shows two simulated electric field distributions with equi-potential lines 
from the sensor with an electrically conductive ring with conductivity ten times as * 
that of the principle flow/material to be'measured; 

Figure 6 shows an extension of the invention using a resistor network for a sensor 
with discrete electrode construction; 
0 Figure 7 shows a sensing strategy with a ti/2 (90°) separation angles for measurement 
and excitation; 
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Figure 8 is the sensitivity distribution from the voltage measurement opposite to the 
current excitation using the %I2 sensing strategy; 
Figure 9 is the overall sensitivity map using the ti/2 sensing strategy; 
Figure 10 is the sensitivity distribution from the voltage measurement opposite to the 
5 current excitation using the adjacent sensing strategy; 

Figure 1 1 is the overall sensitivity map using adjacent sensing strategy; 
Figure 12 is the definition of the sensor geometry; 

Figure 13 shows a comparison of sensitivity of two sensors at the elements across the 
...diameter; 

10 Figure 14 demonstrates the distribution of the illness measure (condition number); 
Figure 1 5 gives the illness curve (a conductivity ratio of 1 stands for conventional 
sensor); 

Figure 16 shows alternative constructions of the conductive sensor; 
Figure 1 7 is a flow chart of SCG reconstruction procedure; 
15 Figure 18 is reconstructed images from simulated data with large conductivity 
distribution difference, (a) conductivity set-up with two conductivity values, 0.11 and 
0.14 mS/cm ; (b) image obtained from the MSBP, (c) image obtained from the 
EIDORS with 6 steps of solution (Tikhonov regularisation parameter lxlO' 3 ), (d, f) 
images obtained from the SCG with 5 steps of solution and 20 iterations taken in the 
20 GCG for solving each inverse function, (e) reconstruction convergence from (d). 

Figure 19 is reconstructed images for investigating imaging distiguishability. (a) 
conductivity set-up with two conductivity valued, 0.055 and 0.11 mS/cm, (b) image 
obtained from the MSBP, (c) image obtained from the EIDORS with 6-steps of 
solution (Tikhonov regularisation parameter 1x1 0~6), (d, f) image obtained from the 
25 SCG with 20 steps of solution and 20 iterations taken in the GCG for solving each 
inverse solution, (e) reconstruction convergence from (d); 

Figure 20 is a schematic layout of the 20 meters How.loop used for the experiment; 
Figure 21 is reconstructed 2D images in respect to typical air cavity formation in the 
flow loop; 

3 0 Figure 22 is a sequence of cross section images of a human hand in 5 . 1 6mS/ cm brine, 
which were reconstructed using STM algorithm with three steps of inversion; 
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Figure 23 shows 3-D hand images interpolated from the 2-D EIT images given in 
Figure 22 and then iso-surfaced with a cutting value 3 .3mS/cm (Spyglass vl .01). 

The electrically conductive ring sensor basically consists of an electrically conductive 
5 ring 1, a number of electrical contacts 2, a target content 3, a number of electrical 
excitation sources 4 and a number of voltage measurement devices 5 [Figure 1]. 

The electrically conductive ring 1 is made of solid substances such as a kind of 
metals or alternatively, other electrically conductive materials such as conductive 
10 rubber or ceramic-metal with a conductivity value much higher than that of the 
content 3 to be measured. The electrically conductive ring acts as an f continuos 
electrode' to distribute the electrical current flow and generate an electric field for 
mapping the impedance distribution of the content 3 due to applying currents from 
the metallic contacts 2. The metallic contacts 2 are embedded into the electrically 
1 5 conductive ring in good electrical contact with the outside wall of the electrically 

conductive ring. The size of these contacts 2 may be very small since a low common 
voltage on the contacts is produced from the high conductive ring 1 and there is no 
effects of electrode-electrolyte interface between these contacts 2 and the electrically 
conductive ring 1 if the ring is made of solid substance. Therefore, the amount of the 
20 contacts to be installed could be much more than the number of discrete electrodes in 
the conventional sensors. A number of electrical excitation sources, such as the 
currents or voltages 4, are applied to one or more electrical contacts 2. Boundary 
voltages or currents are also measured from the electrical contacts 2. A continuous 
electric field distribution presents along the inner surface of the electrically 
25 conductive ring that is in electrical contact with the target content 3 . The target 

content 3 may be a process fluid with electrically continuous phase, such as that in a 
miscible liquid mixing, .or with electrically disconf inuous phases, such as those in 
oil- water-gas flow or froth formation [figure 2]. Since no discrete electrode directly 
contact to the target content 3, the excitation and measurement will not get into 
3 0 saturation under a certain current value when a part of the inner sensing surface loses 
an electrical contact to the principle flow. Figure 3 illustrates its flange configuration. 
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A more homogeneous electric field distribution can be produced if the conductivity 
of me electrically conductive ring 1 is higher than that of the target content 3 . 
Comparisons of the potential distributions generated from a sensor with discrete 
5 electrodes (Figure 4a and 4b) and the electrically conductive ring sensor (Figure 5a 
and 5b) are given, which were simulated using Ansoft Maxwell 2D simulator 
[Version 1-7-66, 1996, Ansoft Corporation]. In the simulation, the conductivity of 
me electrically conductive ring was ten times as that of the target content. The 
potential distributions in Figure 4a and 5a were generated from a current excitation 

10 between an adjacent electrode/contact pair. Figure 4b and 5b are potential 

distributions relating to an opposite electrode/contact pair excitation position. The 
potential distributions in Figure 5a and 5b are more equal-spaced man those given in 
Figure 4a and 4b. It is well known that the field-stream lines can be derived from the 
equal-potential lines based on the orthogonal method. Therefore, the simulation 

15 demonstrates that the electric field intensity distribution generated from the 

electrically conductive ring sensor are more uniform than that generated from a 
discrete electrode sensor if the conductivity of the ring 1 is high than the conductivity 
of the target content 3 . 

20 With the previous sensor construction [Figure 6], a number of discrete electrodes 6 
were intrusively embedded into anon-conductive vessel 7 with good electrical 
. insulation between each other and electrically contact the target content 3. For a fluid 
with discontinuous phases of electrical conductivity, such as a half filled horizontal 
stratified flow or a foam construction with a large bubbling formation, it could not 

25 provide sufficient measurements since some of electrodes have lost their electrical 
contact with the conductive phase of me fluid [Figure 2]. But for a conductive ring 
sensor, the current is always continuous even in situation like figure 2. As an 
extension of the invention, a resistor network 8 is proposed for such sensors with 
discrete electrode structure to perform a completed set of measurements. These 

30 electrodes are connected with the resistor network 8 using electrically conductive 

wires 9. The resistor network 9 plays a similar role as the electrically conductive ring 
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1, which -bypasses and re-distributes currents in case of the current excitation 
electrodes blocked by a discontinuous fluid or bubbles, and produces a more 
homogeneous electric field than that produced fronra conventional discrete electrode 
sensor. Figure 6 shows the extension for the sensor with discrete electrodes. 

The measurement relationship between one and other sensing strategies for a 4- 
eiectrode measurement protocol can be derived as Equation 1 & 2 based on the 
reciprocity theorem [Geselowitz, IEEE Trans. Biomed. Eng. BME-18], which may be ' 
summarised as 'with a 4-electrode measurement protocol, the boundary 
measurements , and the sensitivity .matrix for an alternative sensing strategy can be 
derived from the complete set of independent measurements and the sensitivity 
matrix obtained from a known sensing strategy'. The measurements based on the 
alternative sensing strategy can not generate any more new information than that 
given from the complete set of independent measurements, but could produce a' 
different signal-to-noise ratio and a different sensitivity distribution over the imaging 
domain. A general form of the boundary voltage relationship between an alternative 
4-electrode sensing strategy and adjacent sensing strategy can be expressed as 
Equation 1, which derives mutual impedance or boundary voltages for an alternative • 
4-electrode sensing strategy from those for adjacent electrode pair strategy. 

Z ItJ d MM )= EE Z, M (/ m>ra+I ) (Equation 1) 

(1 < M < N <I < J < number of electrode ) 

Where ZjJJmn) represents the mutual impedance obtained from voltage measured 
between electrode / and / when current presents between electrode Mand K 

The sensitivity matrix of an alternative sensing strategy can also be derived from 
algebraic combination of the sensitivity matrix obtained from a complete set of 
independent measurements because a linear approximation is adopted in the 
calculation of the sensitivity matrix, and the linear relationship between the two sets 
of boundary voltages exist (Equation 1). Derivation to a new sensitivity matrix for an 
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alternative strategy from the sensitivity matrix of adjacent electrode pair strategy is 
given in Equation 2. 

K-l 7-1 

s i,j,k(lM,N) = 2u2j S *\m>* C^,«+i) (Equation 2) 

5 (l^ M <N <I < J< number of electrodes ) 

Where S f J (I UN ) represents the sensitivity coeiBcient at pixel, 1c, obtained from the 
measurement and excitation position where the voltage is measured between 
Electrode I and J when current presents between electrode M and N. 

1° 

Based on Equation 2 5 a sensitivity distribution with improved homogeneity, 
especially at the central area, has been found from a sensing strategy that the angle 
between two electrodes/contacts for measurement and excitation is n/2. This n/2 
spaced sensing strategy is demonstrated in Figure 7 3 which gives a total of 72 
15 independent measurements. Figure 8 shows one of the sensitivity maps relating to the 
measurement position opposite to the current injection. Summing the coefficients 
from the seventy-two normalised sensitivity maps at each pixel, it gives an overall 
sensitivity distribution as Figure 9. To compare the sensitivity homogeneity, two 
related sensitivity maps obtained from adjacent sensing strategy are given in Figure 9 
20 and Figure 10. It can be seen that Figure 8 gives a more homogenous sensitivity band 
across the central area, compared to Figure 10. The overall sensitivity distribution in 
Figure 9 is more homogenous as well, compared to that given in Figure 1 1, which 
demonstrates that a much improved total-sensitivity distribution over the sensing area 
has been produced. The results of the alternative sensing strategy imply that the 
25 imaging sensitivity can be relatively intensified or attenuated at a particulate 

interested area in the whole imaging domain. The sensitivity matrix for such specific 
distribution can be derived from either' the alternative sensing strategy or the 
particulate combination- of the sensitivity matrix derived from a known basic set of 
independent measurements (e.g. adjacent electrode pair sensing strategy). 

30 

The sensor geometry and its conductivity value can be further optimised. 
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It is found that the ill conditioning of the sensitivity matrix (essentially Jacobean 
matrix) is related to the non-homogeneity of the sensitivity field. The bigger the 
condition number of the sensitivity matrix is, the more non-homogeneous the 
sensitivity field is. The smallest sensitivity is related to the central elements with 
5 largest related to periphery elements [Yorkey, et al, IEEE Trans. Biomed. Eng., 
BME 34(1 1), 1987]. Take the condition number as the optimum objective, we can 
optimize the conductive ring sensor against different thickness and conductivity of 
the conductive ring sensor. 

1 0 Figure 12 shows the sensor geometry used during the optimum design. Here, sixteen 
electrical contacts are distributed around the conductive ring, possessing half of the ■ 
ckcumference of the external side of the ring. The ring itself has an internal diameter 
of 5cm. The conductive ring thickness, t, was varied during the simulations. The 
value of t was assumed to be 0.2cm, 0.6cm, 0.8cm, and 1 .0cm. The conductivity of 

1 5 the ring is also a variable ranging from 0. 1 lms/cm to 0.66ms/cm. The sensor was 
assumed filled with water with a conductivity of 0 . 1 lms/cm. 

A comparison of sensitivity distribution across the diameter is shown in figure 13, 
where the current is injected from electrode 1 and 2 with measurements taken from 
20 electrodes 9 and 10. The distribution in the conductive ring sensor shows more 
homogenous than that in traditional sensors. 

It is worth to mention that the important parameters are the ratio of the conductivity 
and thickness of the ring to those of the material filling the sensor, rather than their 
25 absolute values. For example, the numerical results for one case could be obtained 
from another case by appropriate scaling if the same conductivity and thickness ratio 
were used in simulations. 

A "sensitivity matrix for a commonly used adjacent sensing protocol [Brown et al, 
30 Proc. IEE Lit Conf. on Electric and Magnetic field in Med. and Bio, 1985] is 
computed and condition number calculated for each conductivity and thickness 
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value. Figure 14 shows the distribution of the condition number of the sensitivity 
matrix when the thickness and conductivity of the conductive ring sensor varies. 

From figure 14, we could find that within a modest range of conductivity, the illness 
5 measure (the condition number) decreases with the increase of the conductivity and 
the thickness of the conductive ring sensor, as was confirmed by visualisation with 
Maxwell 2D simulations. . 

It is important to note that the illness measure is not a monotonously decreasing 
1 0 function of the thickness and conductivity of the conductive ring sensor. As an 
example, when the thickness of the conductive ring sensor is fixed at 1, and the 
conductivity of the conductive ring sensor is a variable, the condition number curve 
is shown in Figure 15. It is found that the illness measure reaches its minimum at the 
conductivity of 0.66ms/cm. In Figure. 14, apparently, the optimal parameters of the 
1 5 conductive ring should be 0.55ms/cm and 4mm for the conductivity and the thickness 
respectively. 

In addition to the traditional "tube" shape sensor configuration, the sensor can also be 
constructed from a conductive disk 1 1. A hole is drilled in the center of Hie disk and 
20 a layer of metallic co ating serving as electrode contacts 2 is made around the ring. 
[Figure. 16a]. 

The disk shape sensor can also be configured into a flange or a washer from 18 
[Figure 16c]. In order to facilitate the wiring, PCB 1 5 can be made to integrate the 
25 electrical contacts 2, wiring 13 and connector 12 [Figure 16b]. Another advantage of 
this configuration is its flexibility for replacing the conducive ring once it is 
contaminated such as in food or pharmaceutical applications. 

Following the idea above, the conductive sensor can be made in a micro- or nano- 
3 0 scale using micro-machine technology or integrated circuit fabrication technologies. ' 
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The sensor can also be fabricated integrated with its processing circuits on a single 
silicon chip. 

The sensitivity theorem deduced by Geselowitz [TEfiE Trans. BME-18, 1971] and 
5 later refined by Lehr [IEEE Trans. BME-19, 1972] for a two ports system can be 
expressed in discrete domain as Equation 3 [Murai, et al, IEEE Trans. BME-32, 
1985] and Equation 4 [Wang, etal, Chem. Eng. Comm. VI 75, 1999% assuming the 
conductivity change is much smaller than the original value (Aoj c « <j$ and 

-neglectjng,theilLigh_ o rder terms and supposing that the conductivity distribution is 
10 composed from w small uniform patches' or pixels. In this case, the sensitivity 
coefficient, s, for each discrete pixel is given by Equation 5 [Murai, et al, IEEE 
Trans. BME-32, 1985], where 77^ is a discrete area at k. 



25 



15 s t 



AZ = jAcr^^ (Equation 3) 

Z = ~YjCr k s^ k (Equation 4) 

V<f> Vyr 



W = ~I ~r F" ^* (Equations) 



where y/ AB , <f> CD are potentials measured from ports A-B and C-D in response to the 
present of currents 1 Y and 1^ respectively. 

In practical applications of electrical impedance tomography, image reconstruction 
20 was highly affected by the electrode modelling error and the measurement noise. In 
practice, it was hard to make the electrode system fully consistent (e.g. geometry, 
electrode-electrolyte interface, installation environment, etc.). Quality of data 
acquired under industrial environment was also decayed by industrial and 
instrumental noises. 



Normalisation procedure has been employed in almost all back-projection 
algorithms. The normalisation procedure, as given in Equation 6, can be derived from 
dividing Equation 3 by Equation 4, 
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C(A* k <<tii,J°*l,2,...>F) (6) 

Wherey is the measurement-projection location and £is the pixel number, s jk denotes 
the sensitivity coefficient at pixel k under the measurement-projection j\ P denotes 
the maximum number of measurements, w denotes the maximum number of pixels, 
- -0-^-and ^^-are4fae^conducthdte_an d conductivity change at pixe l k, respectively, and 
Vj and AV. are in respect to the reference voltage and the voltage change at 
measurement-proj ection j. 

The problem of conductivity distribution is of a nonlinear type. As the sensitivity 
theorem is based on a linear approximation with a condition Aa k «cr k , therefore, an 
iterative approach was employed in many cases to approach the true value. Both 
forward and inverse solution procedures have to be employed in the multi-step 
solution, which are used for error estimation and updating, and resulting a step 
solution respectively. Conjugate gradients (CG) method is utilised to solution the 
linear equations in both forward and inverse problems. In this patent this method is 
named as sensitivity theorem based inverse solution using conjugate gradient 
methods (SCG). 

Forward solution in the multi-step inverse solution is used for producing an error 
vector for each step inversion. It also up-dates the sensitivity matrix for next step 
inverse solution. Sensitivity matrix can be derived directly from the nodal voltages 
obtained from the forward solution [Yorkey, et al, IEEE Trans. BME-34, 1987, . 
Murai, et al, IEEE Trans. BME-32, 1985]. The actual current value used in the 
solution is not significant as long as the value doesn't vary in whole process since 
only relative changes of the boundary voltage measurements are employed in 
following inverse solution. 



J k=\ 

y ~ ™ 
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The solution of a FEM model presenting a 2D cross-section of a process vessel with 
Neumann boundary conditions in addition to a single Dirichlet condition to avoid 
singularity can be solved by Equation 7 utilising a linear approximation given by 
Yorkey et al and Murai et al.[Yorkey, et al, IEEE Trans. BME-34, 1987, Murai, et al, 
5 IEEE Trans. BME-32, 1985]. 

Yv = c (7) 
where Y, v, c denotes the global conductance matrix, the nodal voltage vector, and 
. the nodal current vector resp ectively. 
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Based on the expression given by Eq.6 and the assumption of a homogenous 
conductivity distribution, a 0> at the time of taking reference V, the inverse solution in 
the multi-step approach is given by Eq.8 in the form of matrix notation, 

T - -s ' e (8) 
where the elements of the normalised sensitivity matrix s^ 1 at the iteration n, the 
vector of conductivity relative change ^at pixel k, and the vector of boundary relative 
change e at projection j are denoted as Eq.9, Eq.10 and Bq.ll, respectively. The 
conductivity is updated by Eq. 12. 

J(B) _ ; *jM n) ) 

s /jc ~— r (9) 

2X*(^ 0) ) 



jo* - A^' J 



n cJ?> (10) 



M AV _ r](a)-Vj(<r) 



e y=ir- VM) (ii) 



^) = ^)(1 + ^)) C12) 
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where a and & are actual conductivity distribution at the moment of acquiring the 
reference voltage V and measurement- V. <r< 0) and are estimated conductivity 
values for simulating cj and c/, 

5 Since conductivity is inversely related to voltage, the conductivity updating also can 
be based on an approximation of the inverse relation (Eq.13 based on 1+x & l/(l-x) at 
x<l) 3 which can improve the convergence speed for both positive and negative * 
changes in conductivity. 

10 ^rS 1 w 

Noting the validating condition Acr«cr 3 the s M (<x^) in response to Uj(df )y ) as 

well as the regularisation procedure in the linear approximation Eq.6, Eq.ll thus 
15 decomposes to 

20 

where the and W) are the estimated reference voltage and updated 

measurement voltage in respect to the conductivity distribution <?* 0) and o 1 ** 0 . After 
a number of steps of updating the conductivity distribution and the sensitivity matrix, 
the decomposed boundary relative change or error vector in Eq.14 will be minimised. 
25 It is thought the non-linear inverse solution has been reached when the norm of the 
error vector is sufficiently small. 

The procedure of multi-step iteration is performed as follows. 

30 (1) Pre-compute the assumed boundaiy voltage vector and the sensitivity matrix, 

u(^ } ) 3 s(^) 

(2) Measure the boundary voltage profiles and produce the relative change vector, 
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(3) Pre-set the iteration control factors for the minimum convergence error, e s > the 
maximum number of inverse steps, 5 S , and the maximum number of iterations, S c> 

5 for the GCG. 

(4) Initiate the first estimations for the error function vector and the conductivity 

vector, 

a« =[ £ r«,...,cr«] r 5 . tff =o- 0 /ci:^/ p ) 

7=1 

10 (5) Normalise the sensitivity matrix (the iteration steps n is from 7 to 8 S ) 



E^ 0) ) 
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(6) Solve the inverse problem using the GCG method and then update the relative 
change in the conductivity vector (Maximum Iteration = S c ), 

y = —s ~ l • e 
0 (7) Update the conductivity vector, 

(8) Solve the forward problem using the ordinary CG method to update the boundary 

voltage vector and sensitivity matrix, 

u(c^ +l) ); s(c^ +1) ) 

(9) Update the error vector, 



(1 0_) Check the control factors to determine whether one of them has been reached, ' 

||e)|<5-, or Steps ><5, ? 
17 
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(11) If the convergence or the number of maximum iterations is not reached, the . 
process jumps to step (5) until one of these conditions is reached. To summarise the 
5 procedure, its flowchart is given as Figure 17. 

The sensitivity coefficient back-projection approximation method (SBP), utilising a 
normalised transpose matrix of the sensitivity matrix obtained from Eq.5 as a 
weighting matrix, was first refined for ERT by Kotre [Kotre, Physiology 
10 Measureraentsr45y4fappL^ A-199Zf >_Ihe-ap proximation has been successfully 
applied to many industrial applications [Williams et al, Process tomography — 
principles, techniques and applications,, Butterworth-Heinemann Ltd., 1995] due to 
its good anti-noise capability and fast solution speed. The physical definition may 
follow the basic principle of linear back projection: the relative change of the 
1 5 boundary voltage measurements is linearly back-projected to each pixel between two 
equipotential lines in the case of the equipotential back projection [Baber et al, 
Applied potential tomography', J. Phys. E: Sci. lustrum., 17,1984], or over whole 
domain in the case of the SBP [Kotre, Physiology Measurements, 15, Suppl. 2A 
1992]. The sum of the products of the relative change and the weight 
20 factor/sensitivity coefficient at each pixel, normalised by the sum of their weight 
factors/sensitivity coefficients that are derived from all possible boundary 
measurements, approximately presents the relative change of the conductivity at this 
pixel. 

25 To investigate the accuracy and the limitation of the multi-step solution, a sequence 
of images were reconstructed from simulated data. A part of these results are 
presented in this patent. Only the adjacent sensing strategy [Brown et al, Proc. IEE 
hit. Conf. on Electric and Magnetic field in Med. And Bio., 1985] was employed in 
these simulations. All references were taken from a homogeneous set-up with a 

3 0 conductivity of 0 . 1 lmS/cm. Meshes with 1 04 and 224 triangle pixels were used for 
simulating boundary voltages from these conductivity set-ups. To investigate the 
discretization error and mesh adaptability of the algorithm, some images were 

18 • 



WO 02/053029 



PCT/GB01/05636 



reconstructed using the mesh with 224 elements for data simulated using the mesh 
with 104 element. Electrode positions at the 224 element mesh also has an 1 1.25° 
anticlockwise rotation compared to those in the simulated set-ups with 104 element 
mesh. 

5 

Considering the importance of the SBP in current applications, these simulated data 
were also reconstructed using the SBP. A mesh with 316 square pixels was employed 
for the SBP algorithm. The electrode position at the square mesh has an 1 1 .25° 
anticlockwise rotati oncopipared to those in the simulated set-ups with the 1 04 
1 0 element mesh. 

To assess the imaging capability of the SCG algorithm, a Newton-Raphson method 
based algorithm with Tikhonov regularisation [Vaukonen, Scientific Abstracts of 
2nd EPSRC Engineering Network Meeting, UCL, London, 2000], named as EIDORS 

15 (EIT and Diffuse Optical Tomography Reconstruction Software, UMIST 2000), is 
introduced to reconstruct a part of these data. These images are reconstructed using a 
mesh with 492 elements and parameters of the simulated contact 
impedance=0.005Q/cm, Tikhonov regularizationparameter=0.001 and the iteration ' 
number=6. The electrode position also has a clockwise rotation compared to the 

20 simulation set-up. The original simulated data are directly used for reconstructing 
images with the EIDORS. However, the data obtained from real measurements have 
to be regularised by applying the products of relative changes and the standard 
boundary voltages from a homogeneous set-up in the EIDORS. 

25 A set-up with a complex conductivity distribution, as shown in Fig. 1 8 a, is 
reconstructed as Figs. 1 8b and 1 8c from the SBP and EIDORS algorithms 
respectively. It is obvious that the SBP algorithm could not deliver an accurate image 
for this set-up (Fig. 18b). Both the EIDORS (Fig. 18c) and SCG (Fig.l8d & f) can 
reconstruct the complexity of the set-up. The image obtained from five steps of the 
0 SCG solution (Fig. 1 8c) provides a conductivity error of 3 .32% and a boundary 
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voltage error of 0.055%. The discretisation error caused by applying different mesh is 
reflected in Fig.3f. 

Imaging distinguishability of these algorithms was investigated and the results are 
5 given in Fig. 19. Two 'objects' were located along a radius of a mesh (Fig. 19a). For 
the set-up, SBP can image the presence of the two objects but could not distinguish 
them (Fig.l9b). The edge of the objects' image is also merged with the boundary. 
The image obtained from EDDORS (Fig. 19c) gives a better presentation about the 
object location with a clear separation from the boundary but is incapable of 
0 distinguishing between the two objects. SCG can clearly distinguish between the two 
objects although a certain artificial background noise presents (Fig. 19d & f). 

increasing the number of steps to solve the non-linear problem may improve the 
imaging accuracy. However, the most significant contribution is from the first step, 
> which has already been shown in previous convergence graphs. Therefore, in 
practice, a single step solution would be widely acceptable for the balance of the 
imaging accuracy and speed. 

The single step method solution is be further simplified by obtaining an 
approximation of the inverse of s (s is the normalised sensitivity matrix) and 
therefore the whole computation for reconstructing a single image only involves a 
multiplication of a matrix and a vector. The detailed information is provided as 
following. 

The inverse solution can be denoted by equation (15) in a matrix form 

e = -s • y (Equation 15) 

Where s is the normalized sensitivity matrix, y the conductivity relative change 
vector, and e the boundary relative change vector. 
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Since s is neither square nor positively definite, direct solution to equation (15) 
using CG method would not exist. In order to be solved, equation (15) has to be 
transformed to equation (16) by applying minimization to equation (15). 

5 s r e - -s r sy (Equation 16) 

With s T s presenting a very large condition number, equation (1 6) cannot be solved 
using direct inverse method aiming to an accurate solution. An iterative method with 
-fimte-steps-gives-equation^ which serves as a kind of 

10 regularization. Although being able to provide a sensible solution to equation (16), 
iterative methods involves much larger computation load than direct inverse methods 
by which an inverse matrix for s r s can be obtained, and the solution to equation 
(1 6) only involves a multiplication of a matrix (s r s)" 1 • s T with an vector e . 

15 A.method of obtaining (s T s)~ l using CG method and therefore providing a one-step 
solution to equation (16) is presented, having the advantages of both lower 
computation load and a sensible solution. 

To simplify notations, let A= (s r s)" 1 , and form n equations in (1 7). 





T 




"0" 




"0" 




"0" 






0 




1 




0 




0 




A*xl= 




, A*x2= 




, A*x3= 


1 






. (Equation 17) 




_0_ 




_0_ 




_0_ 




1_ 





Solving equations (17) using CG method with finite steps (typically 20), we can get 

vectors, xl, x2, x3 xn, which can be assembled into [xl, x2, xn], an 

25 approximation of (s T s)~ 1 . 
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Knowing ?, (s T s)" 1 -s T can be pre-calculated and therefore equation (16) can be 
solved in one step. Let T*-(s T J) -1 • s T , then the reconstruction process can be 
denoted as: 

cr«flr 0 *a + r) = ^ 0 *a-r*tf) (18) 

5 

Application 1 

Multiphase flow exists in many industrial processes. The phase distribution and 
-interfaces Jba.ajaiultiph.ase flow carries si gnificant information about the pro cesses 
10 [Chhabra et al, Non-newtonian flow in the process industries - fundamentals and 
engineering applications, Butterworth-Heinemann, Oxford, 1999]. However its 
behaviour is extremely complex, which presents a great challenge to the study of the 
flow mechanisms and the measurement of multiphase flow. Both single modality and 
multi-modality imaging have been reported [Dyakowsfci et al, Powder Techno! 104, 
15 1999; Hoyle et al, Proa 1st World Congress on Industrial Process Tomography, 

Buxton, UK, 1999] as a means to identify flow regimes in two-phase flow. A direct 
imaging method for characterising gas-water flow based on a simple 4-electrode 
technique was proposed [Seleghim et al, Meas. Sci. Technol. 9 1998]. Solid 
concentration distribution can be imaged using electrical resistance tomography 
20 (ERT) and its local velocity profile could be interpreted by applying cross-correlation 
method for a twin-plane sensing system in a full-filled pipe flow [Lucas et al, Journal 
of Flow Measurement and Instrumentation, 10, 4,1999]. However, most previous 
flow applications using conventional electrical resistance tomography were only 
applied for vertical pipes or horizontal pipes in the absence of large gas bubbles. 
Signal saturation will occur for a half-filled pipe flow, such as churn flow in a 
horizontal flow loop, since part of electrodes lose their electrical contacts with the 
electrically conductive fluid. Imaging methods based on the electrical conductive ring 
sensing strategy has been used for measuring the two-phase flow with 'electrically 
continuous' and 'discontinuous' phases. As a preliminary test, several examples of 
identifying flow regimes using this technique are presented 
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A.20 m flow loop (inner diameter, 50mm) fitted with ERT sensors was employed for 
the experiments. A schematic diagram and overview of the flow loop are given in 
Figure 20. Mains tap-water (conductivity = 0.304mS/cm) was used as the liquid 
phase and air, were introduced into the flow loop from a mixing jet-pump. 
5 Measurements were performed at ambient temperature. By controlling the air flow- 
rate at the air inlet of the jet-pump, different flow regimes were generated in the low 
part of the flow loop. 

The experiments wereperformed under different air flow-rates of 0.5m 3 /h 5 - 1 .5m 3 /h, 
10 2.5m 3 /h, 3.5m 3 /h and the one less than 0.5m 3 /h in regard to the productions of slug 
flow, slug-churn flow, churn flow and bubbly flow regimes. The water flows were 
scaled with an accumulating tank during the experiments to get water flow rate and 
mean velocity. At the mean time, a number of photographs were recorded as visual 
presentations of these different flow regimes. 

15 

Two-dimensional (2D) -slice images were sequentially reconstructed. Some 
- reconstructed images in respect to typical air cavities in the flow loop are shown in 
Figure 21, where the blue, areas represent the air cavities or low conductivity regions 
and the red areas represent the water or high conductivity regions. 

20 

Application 2 

As an interesting test, a human hand was scanned using an EIT system. A 14.8cm 
diameter vessel fitted with one 16-electrode ring sensor and filled with 5.16mS/cm 

25 brine was employed. The adjacent electrode sensing strategy was applied with a pair 
of 10mA@9.6kHz sinusoidal currents for the test. The image data was made up of 15 
measurement data sets acquired as the hand moved along the axial direction of the 
phantom while the water volume was kept the same. The signal-to-noise-ratio (SNR, 
here is the repeatability of reference voltage) was checked, which was within a 

30 maximum error 0.5%. Fifteen 2D images were reconstructed with the SCG [Figure 
22]. These fifteen 2D SCG images were further interpolated to a 3D hand image. A 
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cutting value of 3.3mS/cm was applied to extract the surface of the hand skin as an 
iso-surface using Spyglass vl .00. An outline of the human hand has been 
successfully extracted as shown in Figure 23. Some distortions can also be found, 
which may be caused by the electrode noise, the unstable position of the hand during 
5 the scanning and the 3-D effect of electrical field. 
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Claims 

1 . A method of measuring the distribution of electrical impedance of a 
multi-phase flow, of which the principal flow has an electrically 

5 continuous or discontinuous phase, using an electrically conductive ring 

electrode, characterised in that the ring electrode uses a 7c/2 sensing 
strategy. 

2. _A method according to Claim 1 ch aracterised in that the ring electrode has 
0 a much higher conductivity compared to the principle flow or materials to 

be measured and it is employed as a ring electrode of a part comprised as 
part of a sensor body. 

3 . A method according to Claim 1 charaterised in that a more homogeneous 
5 electric field distribution is produced by maintaining the conductivity of 

the electrical conductive ring electrode at a higher level than that of the 
multiphase flow. 

4. A method according to Claim 1 characterised in that the electrically 
conductive ring is made of solid substances selected from metal, 
conductive rubber and ceramics or a combination of the aforesaid. 

5. A method according to Claim 1 characterised in that a number of 
electrical contacts are embedded into the conductive ring electrode in 
good electrical contact with an outside wall of the ring electrode or 
penetrated through the wall of the ring electrode in good electrical contact 
with the wall. 

6. A method according to Claim 1 characterised in that the size of the 
contacts is small. 
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7. A method according to Claim 1 5 characterised in that the sensitivity 
distribution of the voltage measurement is relatively intensified at a 
specific area in the whole of the sensing domain. 

5 8, A method according to Claim 1 characterised in that the intensified 

sensitivity matrix is derived from the specific sensing strategy, or a 
combination from a basic sensitivity matrix. 

9. A method according to claim 8 characterised in that the basic sensitivity 
1 0 matrix is derived from equations 1 and 2. 

10. A method according to Claim 1 characterised in that a 7t/2 spaced sensing 
strategy with a 90° separation angle for measurement excitation is 
employed to further improve the homogeneity and the relative intensity of 

1 5 the sensitivity distribution at the central area of the sensing domain. 

11. A method according to Claim 1 characterised in that the target is selected 
from a process material, a process fluid and a part of a human body. 

20 ' 12. A method according to claim 10 characterised in that the electrically 

continuous phase is a miscible liquid or the discontinuous phase such is in 
a foam formation or a half filled horizontal oil/water/gas stratified flow. 

13 . A method according to Claim 1 characterised in that the method is 

25 extended to a sensor with discrete electrode structure using an external 

resistor network connecting all electrodes. 

14. A method according to Claim 1 characterised in that the geometry and the 
conductivity of the conductive ring sensor is optimized using the illness 

3 0 measure of the sensitivity matrix. 
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15. A method according to Claim 1 characterised in that a conductive disk is 
used to construct the sensor 3 with PCB as its connection. 

16. A method according to Claim 1 characterised in that the conductive 
sensor is made in micro-scale. 

17. A method according to Claim 15 characterised in that the micro-scale 
conductive sensor is made on silicon chips. 

18. A method according to Claim 1 6 characterised in that the micro-scale 
conductive sensor is integrated with a measurement and processing 
circuit. 

19. A method according to Claim 1 characterised in that a multi-step image 
reconstruction algorithm using an error processing method is employed to 
reconstruct the conductivity distribution. 

20. A method according to Claim 18 characterised in that equation 8 is used. 

21 . A method according to Claim 1 8 characterised in that the error processing 
method uses the difference between the relative changes of the forward 
solutions and the measured results. 

22. A method according to claim 20 characterised in that equation 14 is used. 

23. A method according to Claim 18 characterised in that the CG method is 
used to solve linear equation in implementing both inverse solution and 
forward solution. 
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24. A method according to Claiml 8 characterised in that an image solution is 
performed using a single step to prove an approximated conductivity * 
distribution. 

5 25. A method according to claim 23 characterised in that equation 1 7 is used. 

26. A method according to Claim 21 characterised in that the single step 
method involves only a multiplication of a matrix and a vector. 

0 27. A method according to claim 25 characterised in that equation IS is used. 

28. An apparatus adapted to measure the distribution of electrical impedance 
of a multi-phase flow of which the principal flow has an electrically 
continuous or discontinuous phase, using an electrically conductive ring 
electrode characterised in that the ring electrode is adapted to use a tc/2 
sensing strategy. 

29. A computer program product adapted for determining the internal 
structure of a body with electrical impedance tomography characterised in 
that the programme includes the use of the algorithms as given by the 
equation 8 and 17. 

30. A method of determining the internal structure of a body which utilises a 
method according to Claim 1. 
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